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Abstract 

A recent analysis by de Barros and Suppes of experimentally realizable GHZ cor- 
relations supports the conclusion that these correlations cannot be explained by 
introducing local hidden variables. We show, nevertheless, that their analysis does 
not exclude local hidden variable models in which the inefficiency in the experi- 
ment is an effect not only of random errors in the detector equipment, but is also 
the manifestation of a pre-set, hidden property of the particles ("prism models"). 
Indeed, we present an explicit prism model for the GHZ scenario; that is, a local 
hidden variable model entirely compatible with recent GHZ experiments. 

Key words: local hidden variable, GHZ experiment, detection efficiency, prism 
model 
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1 Introduction 



De Barros and Suppes [1] give a general analysis of realistic experiments, 
where experimental error reduces the perfect correlations of the ideal GHZ 
case. Their analysis makes use of inequalities which are said to be "both 
necessary and sufficient for the existence of a local hidden variable" for the 
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experimentally realizable GHZ correlations. In applying their analysis to the 
Innsbruck experiment [2], however, they only count events in which all the 
detectors fire. While necessary for the analysis of that experiment, they recog- 
nize that this selective procedure weakens the argument for the non-existence 
of local hidden variables. Here we show that they are right and that their anal- 
ysis does not rule out a whole class of local hidden variable models in which 
the detection inefficiency is not (only) the effect of the random errors in the 
detector equipment, but it is a more fundamental phenomenon, the manifes- 
tation of a predetermined hidden property of the particles. This conception of 
local hidden variables was suggested in Fine's prism model [3] and, arguably, 
goes back to Einstein (See [5] Chapter 4). 

Prism models work well in case of the EPR-Bell experiments. The original 
model applied to the 2x2 spin-correlation experiments and was in complete 
accordance with the known experimental results. There appeared, however, 
a theoretical demand to embed the 2x2 prism models into a large n x n 
prism model reproducing all potential 2x2 sub-experiments. This demand 
was motivated by the idea that the real physical process does not know which 
directions are chosen in an experiment. On the other hand, it seemed that in 
the known prism models of the n x n spin-correlation experiment the efficien- 
cies tended to zero, if n — > oo, which contradicts what we expect of actual 
experiments. This contradiction was recently resolved in [10] and [11], which 
show that there is a wide class of physically plausible oo x oo prism models 
with high efficiency (< 0.82). 

In the first part of this paper we explain the principle difference between the 
prism models and the local hidden variable models to which de Barros and 
Suppes' analysis applies. In the second part, we present an explicit prism 
model for the GHZ scenario, a local hidden variable model that is entirely 
compatible with recent GHZ experiments. 



2 The GHZ experiment 



Greenberger, Home, Shimony and Zeilinger [7] developed a proof of the Bell 
theorem without using inequalities. For the GHZ example consider three en- 
tangled photons flying apart along three different straight lines in the horizon- 
tal plane (Fig. 1). Assume that the (polarization part of the) quantum state 
of the three-photon system is 

* = ±= (\H)i ® \ H h ® \V) 3 + l*0i ® \V) 2 ® \H} 3 ) (1) 
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Fig. 1. A three-particle beam-entanglement interferometer 

One can transform the polarization degree of freedom into the momentum 
degree of freedom by means of polarizing beam splitters (see [12]). So the 
quantum state of the system can be written also in the following form: 
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\a) 1 <g> \b) 2 <g> |c) 3 + \a') 1 <g> \b') 2 <g> |c') 3 ) 



where |a) x denotes the particle 1 in beam a, etc. A straightforward interfer- 
ometric calculation ([7], p. 1141) shows that the probabilities of detections 



are 



= 1(1 + sin (0i + 02 + 3 )) 

P2-B+C+ W>1> 02,0a) = |(1 -Sin (0i + 02 + 03)) 

etc. 



PA+B+C+ m,<P2,<P3 



(2) 



(If the number of minuses on the detector labels is even, there is a plus sign; 
if odd, there is a minus sign.) Introduce the following result functions 



1 if the detector A + fires 
— 1 if the detector A~ fires 



B (02) and C (03) have the same meaning for particles 2 and 3. One can also 
show that in state \1/ the expectation value of the product of the three outcomes 
is 



E {A (0i) B (0 2 ) C (0 3 )) = sin (0i + 02 + 3 ) 
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Consider the following choices of angles: 



fti = A (7r/2) B (0) C (0) 
fi 2 = A (0) B (vr/2) C (0) 
fi 3 = A (0) B (0) C (tt/2) 
fi 4 = A (tt/2) B (tt/2) C (tt/2) 

In this case we obtain perfect correlations: 

E(n 1 ) = E(n 2 ) = E(n 3 ) = i (4) 
E(n 4 ) = -i (5) 

So far this is standard quantum mechanics. One can make a Kochen-Specker/EPR- 
type argument, however, if one assumes that in \1/ predetermined values, re- 
vealed by measurement, are assigned to the six observables 

A (tt/2) , A (0) , B (tt/2) , B (0) , C (tt/2) , C (0) (6) 

By virtue of (4) and (5) these values have to satisfy the following constraints: 

Qi = A (tt/2) B (0) C (0) =1 
fi 2 = A (0) B (tt/2) C (0) =1 
fi 3 = A (0) B (0) C (tt/2) =1 
fi 4 = A (tt/2) B (tt/2) C (tt/2) = -1 



Then a contradiction is immediate if we take the product of equations (7). 
Each value appears twice so, whatever the assigned values are, the left hand 
side is a positive number, whereas the right side is —1. 



3 De Barros and Suppes' inequalities 



De Barros and Suppes approach the above contradiction in the following way. 
Without loss of generality, the space of hidden variable can be identified with 
O = {+, — } 6 , the set of the 2 6 = 64 different 6-tuples of possible combinations 
of the values (6). Then the GHZ contradiction amounts to the assertion that 
no probability measure over O reproduces the expectation values (4) and 
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(5). De Barros and Suppes demonstrate this by concentrating on the product 
observables (Qi, . . . , f^) for which they derive a system of inequalities that 
play the same role for GHZ that the general form of the Bell inequalities 
do for EPR-Bohm type experiments [4]; namely, they provide necessary and 
sufficient conditions for a certain class of local hidden variable models. The 
first of their inequalities is just 

-2 < S(fti) + E(Sh) + E(n 3 ) - E(fl 4 ) < 2 

and clearly this is violated by (4) and (5). Moreover if, due to inefficiencies 
in the detectors or to dark photon detection, the observed correlations were 
reduced by some factor e; that is 

E(n 1 ) = E(Sh)=E(n 3 ) = l-e (8) 
E(n 4 ) = -l+e (9) 

then, it follows immediately from this inequality that, "the observed correla- 
tions are only compatible with a local hidden variable theory" if e > |. De 
Barros and Suppes made a detailed analysis of the detection errors and the 
dark photon detections in the realistic detector equipments, and concluded 
that these phenomena do not yield such a large e. 

As in the case of the Bell inequalities, however, the de Barros and Suppes 
derivation starts with the assumption that the variables in (6) are two val- 
ued (either +1 or —1). Since we consider the detection/emission ratio as of 
more fundamental origin, in the prism models developed in the next sections, 
the variables can take on a third value, "D", corresponding to an inherent 
"no show" or defectiveness. In the Bell-EPR case we know that the existence 
of local hidden variables of this more general type are governed by a differ- 
ent system of inequalities. For the inversion symmetric 2x2 case inequalities 
providing necessary and sufficient conditions for prism models were derived 
in [6]. We do not have a comparable system characterizing prism models for 
GHZ type experiments but we will show that GHZ experiments can be mod- 
eled by just such local hidden variable theories. Indeed we will give an explicit 
prism model for a GHZ experiment (with perfect detectors and with zero dark- 
photon detection probability). We will also show that our model is completely 
compatible with the results measured in the Innsbruck experiment. 



4 A toy prism model of the GHZ experiment 

The prism model of the GHZ experiment is a local, deterministic hidden vari- 
able theory, in which the hidden variables predetermine not only the outcomes 
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of the corresponding measurements, but also predetermine whether or not an 
emitted particle arrives to the detector and becomes detected. Consequently, 
the space A of hidden variables ought to be a subset of {+, — , D} 6 . Each 
element of A is a 6-tuple that corresponds to combinations like 

(A(tt/2), A(0), S(tt/2), 5(0), C(tt/2), C(0)) = (+-£>-++) 



which, for example, stands for the case when particle 1 is predetermined to 
produce the outcome +1 if 4>i = 7r /2, —1 if angle 0i = in the measurement, 
particle 2 is n/ 2- defective, i.e., it gives no outcome if 2 = 7r /2, but produces 
an outcome —1 if 02 — 0, particle 3 produces outcome +1 for both cases. The 
essential feature of this conception of hidden variables is that the "values" 
A\(tt/2), A\(0), B\(ir/2), . . . are "prismed" in the sense that, formally, a new 
"value" is introduced, "D", corresponding to the case when the particle is 
predetermined not to produce an outcome. 

Each GHZ event will be represented as a subset U of A. For instance 

£/{A(7r/2)=+}&{B(0)=-}&{C(0)=-} 



stands for the triple detection A + B C with angles (n/2,0,0). 

We have seen that, if determinate values are assigned to all the observables, 
quantum mechanics yields contradictory correlations (7) among the measure- 
ment outcomes at the three stations. Although these four correlations are 
enough for the stated GHZ contradiction, our hidden variable model must be 
consistent with quantum mechanics in a wider sense: the probability measure 
p we are going to define on A C {+, —,-D} 6 must satisfy further constraints, 
following from (2): 

P (U{A(x)=i}k{B(y)=j}k{C(z)=k}\Ul^%) 

= | [1 + ijk sin(x + y + zj\ 



(10) 

x, y,z = 5,0 
for all 2 



i,3,k = ±1 



where U^ z) = U {A (x)^D}&i{B(y)^D}k{c(z)^D}, the event of triple detection. 

Constraints (7) correspond to the fact that some of these probabilities are zero, 
which rule out a large number of 6-tuples. One can show (and easily verify by 
computer) that from the 3 6 = 729 elements of {+, —,-D} 6 there remains 409 
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which satisfy (7). For example: 

• ( — h H — DD) is allowed, because, in this case, whatever the chosen exper- 
imental setup, there is no detection at station 3, consequently there is no 
triple coincidence detection. 

• (D DD+) is allowed because for any measurement setup either the 

outcome triad satisfies the constraints or there is no triple coincidence at 
all. 

• ( D-\ — ) is not allowed, because if the chosen angles were (vr/2, tt/2, tt/2) 

then the results would be A(ir/2) = — , B(7r/2) = — and C(jr/2) = +, which 
would contradict the constraint Q4 = —1. 

There is a prism model on the hidden variable space consisting of these left 409 
elements. However, in order to achieve better detection/emission efficiencies, 
and also to simplify the model, we will refine A further. The 409 combinations 
form four disjoint subsets: 217 of them correspond to the situation where there 
is no triple detection at all, regardless of the angles chosen at the three stations; 
48 combinations produce a triple detection coincidence at only one, and 96 
at two triads of angles (these 48 and 96 form a prism model for GHZ all by 
themselves) and the remaining 48 combinations produce a triple coincidence 
with four different triads of experimental setups. Clearly we achieve the best 
efficiency if we take for A the fourth subset, listed in Table I, and simply omit 
all the others.^ 

For instance the event U B(0) = +" corresponds to 

U{B(0)=+} — (-^3) A5, As, -^10) A15, A17, A20, A22, A24, 

A26) A28, A31, A36, A38, A40, A43, A45, A48) 

Similarly, the event, for example, that u A(tt/2) = + and -B(O) = +" is repre- 
sented by the following subset: 

U{a(it/2)=+}&{b(o)=+} = (A31, A33, A 36 , A 38 , A 40 , A43, A45, A 48 ) 
while, for instance, U{ A (-k/2)=+}&{b{o)=+}&{c{o)=-} = 0- 

Notice that each subset U{ A ( x )^D}k{B( y )^D}k{C(z)^D} - where x,y,z = tt/2 or 

- consists of exactly 24 elements of A. These subsets correspond to the triple 
measurement events that enter into GHZ. 

The probability measure on A must be defined in such a way that the rest of 
conditions (10) (right hand side is not zero) be satisfied. One can verify by 

1 See [8] and [9] for a different presentation of the 48-model, and also for important 
calculations concerning error bounds for GHZ type experiments. 
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Ai= 


— — 


--D+) 


A25 = 


[D + ) 


A 2 = 


— — 


-D-+) 


A26 = 


(£> + - + +-) 


A 3 = 


v 


- + -D) 


A27 = 


[D + + - -+) 


A 4 = 


^ — — 


D-++) 


A28 = 


[D + + + ++) 


A 5 = 


— — 


D + -) 


A29 = 


(+ +D) 


A 6 = 


— — 


+ -+D) 


A30 = 


{+ - -D + +) 


A 7 = 


— — 


+D + -) 


A31 = 


(+-- + £>+) 


A 8 = 


— — 


+ + D-) 


A32 = 




A 9 = 


[-D 


+) 


A33 = 


[+-D + -+) 


Aio= 


[-D 


- + - -) 


A34 = 




A n = 


[-D 


+ - + +) 


A35 = 


(+-+D--) 


Al2 = 


{-D 


+ + + -) 


A36 = 




Al3 = 


{- + 


D) 


A37 = 




Al4 = 


{- + 


-D - -) 


A38 = 


[+D - + + +) 


Al5 = 


(- + 


- + £>-) 


A39 = 


[+D + ) 


Al6 = 


(- + 


D - -+) 


A40 = 


[+D + + -+) 


Al7= 


(- + 


D + +-) 


A41 = 


(+ + - - D-) 


Al8 = 


(- + 


+ -D+) 


A42 = 


+ + 


Al9 = 


{- + 


+D + +) 


A43 = 


(+ + - + +£>) 


A20 = 


(- + 


+ + +D) 


A44 = 


+ D ) 


A 2 i= 


[D - 


- - ++) 


A45 = 


+ D + 


A22 = 


[D - 


- + -+) 


A46 = 


^+ + + - -D) 


A23 = 


[D - 


+ - +-) 


A47 = 


[+ + +D-+) 


A24= 


[D - 


+ + ) 


A48 = 


[+ + + + D+) 



Table 1 



computer that the uniform distribution on A is a suitable one (each element 
has probability ^) and the probability model (A, p) thus obtained has maximal 
triple detection efficiency Indeed, the triple efficiencies are: 



^(triple coincidence) = p (U { A( x )^D}k{B(y)^D}8z{c(z)^D}) = — = 0.5 
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The only way to increase the efficiency would be to modify the probability dis- 
tribution over A. Assuming, however, that for such a non-uniform distribution 
the triple coincidence efficiency is still independent of the chosen experimental 
setups, we have 



48 

^(triple coincidence) = y^p(triple coincidence|Aj)p(Aj) 

i=l 

48 4 1 

= E^p(A i ) = - 

independently of the actual probability distribution p(Aj). 

The key idea of a prism model now is to retrieve the quantum probabilities q(.) 
as the A space probabilities conditional on the measurement outcomes being 
nondefective. Due to conditions (10) this feature of the model is automatically 
provided. Assume, for example, that the chosen angles are {7r/2,0,0}, then 



q({A(n/2) = -}) 

= p(U{A(n/2)=~}\U {A (n/2)^D}&{B(0)^D}&{C(0)^D}) 
P {U{ A (-n/2)=-} n f/{A(7r/2)^D}&{B(0)^D}&{C(0)^D}) 
P (U{A(n/2)^D}k{B(0)^D}k{C(0)^D} ) 

_ P ({Ai, A4, A5, As, Ag, Aio, An, A12, A15, Ai6, A17, Xis}) 

~ 24 

48 

48. _ _ 

24 9 
48 Z 

Similarly, all the other observed single detection probabilities are |. 
To illustrate that constraints (10) are satisfied, consider, for example, 



q{{A{*/2) = +}k{B(n/2) = +}&{C(0) = -}) 

=p (U{A(n/2)=+} n u {Bin/ 2) =+} n u { c(o)=~}\u^ p i e fi) ) 
_ p (U{A(n/2)=+} n u {B(n/2)=+} n U{c(o)=-} n *^f|%) 

fjjtriple \ 
_P({A 3 4, A35, A39}) _ ^ _ 1 



24 






24 


48 






48 






7T 


»)] 


1 — sin 




+ 2 + 




V2 
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Finally, due to the selections involved in building the hidden variable space 
A the model correctly reproduces the GHZ correlations (7), whenever a triple 
detection coincidence occurs: For example the observed expectation value of 
fix is 



E(n,)= E ^(ama|^ 0) )= E -M e , 

, jjtnple , jjtriple P I U/tt q qn I 

A€U (f,o,o) Aety (f,o,o) v i 2 ' u ' u V 



and, similarly, 



E(fi 2 ) = E(Q 3 ) = 1 

E(n A ) = -i 



According to the key idea of a prism model, the above expectation values are 
calculated on sub-ensembles of the emitted particle triads that produce triple 
detection coincidences. In this respect the prism model mirrors actual GHZ 
experiments. 



5 A complete infinite prism model for the GHZ experiment 



In the derivation of the GHZ contradiction we consider only 2x2x2 different 
experimental setups: At each station one considers two possible phase shift 
angles: 0i = |,0; 02 = f,0; 03 = ^,0. In reality, however, the three angles 
can be chosen arbitrarily, and, according to quantum mechanics, the resulted 
triple detection probabilities are given in (2). Although the particular 2x2x2 
scenario suffice to derive a negative statement, the alleged contradiction be- 
tween the existence of a local hidden variable theory and the observed triple 
detection probabilities, it is not sufficient for a positive statement about the 
existence of such a local hidden variable theory. The reason is that if nature 
works according to such a local hidden variable model, then the model must 
reproduce all observed triple detection probabilities for all possible combina- 
tions of angles 0i, 02, 03, since the underlying physical process does not know 
about the angles chosen by the laboratory assistants at the three stations. 
That is why we call a toy model the one we constructed in the previous sec- 
tion, and the existence of a complete infinite prism model covering the whole 
GHZ scenario is still an open question. 

Now we are going to construct a complete infinite prism model for the GHZ 
experiment, which covers the continuum case, when the three phase shift an- 
gles can take arbitrary values. For the sake of later convenience introduce the 



10 



following new parametrization of the phase shift angles: 



a = 0i - g 

7 = 03 - f 

Let us sum up what must be represented in the model: 

I. A continuum set of events corresponding to the detection events, 

together with the triple conjunctions 

At A B+ A C+ 

A+ABtA C- 

P 7 o,U,76 [0,2tt] 

A" A fl£ A C" 
And the "non-defectiveness" events 
A a = A^ V A a 

Bp = B+V B p a, [3, 7 G [0, 2vr] 
C 7 = C+ V C~ 

together with the algebraic relations following from (11). 
II. The "quantum" probabilities of the above events, 

p(^)j(4) jM 

p(A a ) p(A a ) p{Bp) ' p(B ) 
V v (c-) 1 



~ p(C 7 ) ~ p(C 7 ) 2 

p(A+/\B+/\C+) _p(A-ABpAC+) 
p (A a A Bp A C 7 ) p (A a A Bp A C 7 ) 
p(A-AB+AC+) _p(AtABpAC- 



p (A a A Bp A C 7 ) p (A a A Bp A C 7 ) 
= 1(1 -cos (a + /3 + 7 )) 



Fig. 2. The hidden variable space is the union of eight regions G +++ , G ++ ,...G 
The first one is shown in the figure 

p(A-AB+AC+) _ p(A+M3p/\C+) 

p (A a A Bp A C 7 ) p (A a A Bp A C 7 ) 
p{A+AB+AC-) _ p(A-ABp AC-) 

P (A* A Bp A C 7 ) ~ p(A Q A^AC 7 ) 

= i(l + cos(a + /5 + 7 )) (14) 

in accordance with (2). Probabilities p (A a ), p (Bp),p (C 7 ) andp (A a A Bp A C 7 ) 
are the "single and triple detection efficiencies" . 
III. The obvious symmetry of the whole experimental setup: None of the three 
stations is privileged. That is, 

p (Aa) = p (Bp) = p (C 7 ) = uj = constant (15) 



p (A a A Bp A C 7 ) = p (Ap A B a A C 7 ) 

= p(A 1 ABpAC a ) = ... (16) 

Now, the local hidden variable model we are going to construct is based on the 
hidden variable space consisting of eight regions G~ ++ , G + ~ + G . 

Each one is a space S 1 x S* 1 x S 1 , represented by a cube of size (2n) x (2ir) x (2n), 
in which the points of coordinate and 2ir are identified. (The first such region 
G +++ is shown in Fig. 2.) The normalized probability measure is given by the 
eight non-negative density functions p +++ , ■ ■ ■ p , such that 

2tt 2tt 2tt 

J J J p +++ (x,y, z)dxdydz 



2tt 2tt 2tt 

+ J J J p~ ++ (x,y, z)dxdydz 
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2tt 2tt 2tt 

...+ J J j p (x,y, z)dxdydz = l 
ooo 

The events are represented in the following way: 



A + a = { 




G | 


a^rr^a + Aj- 


U \ 


(x,y,z) 


G | 


a<x<a + A}- 


U \ 




G G + — | 


a<x<«+A> 


u \ 


(x, y, z) 


G G + ~ | 


a < x < a + A } 


K = \ 


(x,y, z) 


G | 


a<a;<a + Aj> 


u \ 


(x,y,z) 


G G- + - | 


«<x<a+A| 


u j 


(x,y,z) 


G G~ + | 


a<x<a; + A}- 


u \ 


(x,y,z) 


G G- - | 


a < x < a + A\ 




(x,y,z) 


G C++- | 


7 < z < 7 + A 1 


u j 


(x,y,z) 


G | 


7<z<7+A| 


u \ 


(x,y,z) 


G G+- | 


7<2<7+Al 


u j 


(x,y,z) 


G G | 


7< z < 7 + A > 








a; < x < a + A 


A+A J B+AC+=< 


(x,y,z)eG +++ 


/9<l/</9 + A > 








7 < z < 7 + A 








\ 

a < x < a + A 


^A5+AC+=< 


(x,y,z) EG ++ 


/9<y</9 + A » 








7 < z < 7 + A 


A a = { 


(x,y,z) 


G G+++ | 


«<x<a+A| 




(x,y,z) 


G G- ++ | 


a<a;<a: + A} 




(x,y,z) G G | 


a<x<a; + A| 


C 7 = { 


(x, y, z) G G +++ |7<^<7 + a} 
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U {(x, y, z) G G~ ++ |7<2<7 + A} 

U {(x, y, z) G G |7<2<7 + A} 

One can easily verify that the required symmetries (15) and (16) are satisfied 
by the following Ansatz: 

p +++ (x, y, z) = p~~ + (x, y, z) = p~ + ~(x, y, z) 

= p + ~~(x, y, z) = f(x + y + z)p(x + y + z) (17) 
p~ ++ (x, y, z) = p + ~ + (x, y, z) = p ++ ~(x, y, z) 

= P~~~(x, V, z)= (^--f(x + y + z)^j p(x + y + z) (18) 

where p and / are arbitrary non-negative functions satisfying the following 
conditions: 



2tt+A 2tt+A 2tt+A 



J J J P( x + V + z)dxdydz = 1 (19) 



0</H<^ we[0,6n] (20) 



Now, the probability measure on the hidden variable space, that is, the func- 
tions p and / must be defined in such a way that the quantum probabilities 
(12)-(14) are reproduced. Due to (17)-(20), equation (12) is automatically sat- 
isfied, and if p and / satisfy (13) then they automatically satisfy (14). So, there 
remains only one equation to be solved: 



J 7 7+A C A Ia +A f(* + y + fM; + v + z)dxd y dz 

jjr+A jfj+A ja + A p{x + y + z)dxdydz 

= i(l-cos(a + /3 + 7 )) (21) 

So we are looking for non- negative real functions p(w) and f(w) defined on 
the interval [0, 6ir], satisfying (21) and the conditions (19) and (20). We know 
that f(w) = if w G Ufe=o i 2 3 [2A;7r, 2/c7r + 3A] because cos (2kn) — 1, k — 
0,1,2,3, and f(z) = \ if w G Ufc=o,i,2 [(2A; + l)vr, (2A; + 1)tt + 3A] because 
cos((2/c+ 1)tt) = — 1, k = 0,1,2. These two regions must be disjoint, con- 
sequently A < |, which means - in this model - a limitation for the single 
detection efficiency: u < |. Let us chose A = 0.9|, uo = 15%. 

One can solve (21) numerically. Figure 3 and 4 show the numerical solutions 
for functions / and p. Figure 5 illustrates the high precision of the numerical 
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Fig. 3. A numerical solution of the integral equation (21) for function / 



0.002 






fv 




S 10 12 14 16 IB 



Fig. 4. A numerical solution of the integral equation (21) for function p 




Fig. 5. This is two curves in coincidence: One is the numeric result of the left hand 
side of equation (21), the other is the function on the right hand side 




Fig. 6. The dependence of the triple detection efficiency on a + /3+j. The horizontal 
line represents the value = 0.34% corresponding to the case of independence 

solution. The triple detection efficiency depends on the phase shift angles. 
Figure 6 shows the dependence of the triple detection efficiency on a + f3 + 7. 
The minimal triple efficiency - in this example - is about 0.2%. 
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6 Recent experiments 



Figure 7 shows the schematic drawing of the experimental setup of the Inns- 
bruck experiment [2]. With a small probability, an UV pulse causes a double 
pair creation in the non-linear crystal (BBO). The two pairs created within the 
window of observation are indistinguishable. It can be shown that by restrict- 
ing the ensemble to the sub-ensemble of cases when all of the four detectors, 
T, Di,D 2 , D 3 fire, we obtain the following quantum state: 

i= (\H)i ® \H) 2 ® \V) 3 + \V), ® \V) 2 ® \H) 3 ) ® \H) T 



where \H) T denotes the state of the photon at detector T. This quantum state 
corresponds to a four-particle system consisting of an entangled three-photon 
system in GHZ state, and a fourth independent photon. So we may assume 
that the statistics observed on the sub-ensemble conditioned by the four-fold 
coincidences are the same as those taken on the sub-ensemble conditioned by 
the triple detections at D X ,D 2 and D 3 . What is important from our point of 
view is that any further experimental observations testing the GHZ correla- 
tions, which are based on the above described preparation of GHZ entangled 
states, will be performed on selected sub-ensembles conditioned by the triple 
coincidence detections. Therefore, all of these experimental observations will 
be treated by our local hidden variable model. 




Fig. 7. The experimental setup for demonstration of GHZ entanglement for spatially 
separated photons 

Finally, notice that triple detections, where permitted by the prism model, are 
subject to ordinary sorts of external detection error. If the external detection 
efficiency is, say, d, then triple outcomes having probability 

^(triple detection | none are defective) = 1 



according to the ideal case specified in the model, will have a reduced probabil- 
ity of d 3 , as in the usual analysis of random errors. Similarly we can take into 
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account the non-zero probability of random dark photon detections and make 
a calculation like that of de Barros and Suppes, resulting in the modified ex- 
pectation values (8) and (9). Thus our local hidden variable framework allows 
for the usual techniques of error analysis to treat experimental inefficiencies 
reflected in the actual observations. 

The research was partly supported by the OTKA Foundation, No. T025841 
and No. T032771 (L. E. Szabo). 
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